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Heavy quarkonium radiative transitions have been studied for long time 
in the framework of phenomenological models 114 (for a recent review 
see 15 ). In a typical set up , the starting point is a relativistic Hamiltonian 
whose dynamical fields are heavy quarks and photons; gluons are integrated 
out and replaced by a scalar and a vector interaction. The non-relativistic 
nature of heavy quarkonia is exploited by performing a Foldy-Wouthuysen 
transformation. Anomalous magnetic moments are added as if the system 
was made of two particles independently interacting with the external elec- 
tromagnetic field. 

Such an approach, despite phenomenological success, has an obvious 
limitation: its connection with QCD is unclear. More precisely, while it 
relies on a systematic (non-relativistic) expansion to describe the electro- 
magnetic interaction of the heavy quarks, it does not to describe their 
strong interaction. As a result, the nature of the binding potential between 
the two heavy quarks and of the couplings of the heavy quarks with the 
electromagnetic field remains elusive. In particular, one may wonder if the 
heavy quarkonium anomalous magnetic moment gets contributions from 
the interaction of the heavy quarks in the bound state and if, in the non- 
relativistic expansion, the many couplings between the heavy quarks and 
the electromagnetic field, allowed by the symmetries of QCD, really orga- 
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nize as if they were generated by only two relativistic form factors (i.e. the 
scalar and vector interactions). Finally, the lack of a systematic expansion 
reflects in the fact that model-dependent determinations arc affected by 
unknown uncertainties. 

In the modern framework of non-relativistic effective field theories 
(EFTs) 16-20 , one takes advantage of the small heavy quark velocity v to 
organize calculations of heavy quarkonium observables so that theoretical 
uncertainties can be evaluated and systematically reduced. Applied to ra- 
diative transitions, the EFT approach provides eventually definite answers 
to the above questions 21 . This is done through the following step. 

First, one has to identify the relevant scales in the system. These are 
the heavy quark mass to, the typical momentum transfer in the bound 
state, which is of order mv, and the typical binding energy, which is of 
order mv 2 . The momentum of the emitted photon, fc 7 , is of order tow 2 for 
transitions between different radial levels (like magnctic-dipole hindered 
transitions) and usually of order mv 4 for transitions between levels with the 
same principal quantum number (like magnetic-dipole allowed transitions). 
Since the typical distance between the two heavy quarks is r ~ 1/ (mv) , then 
fc 7 r < 1 and the external electromagnetic field can be multipole expanded 
(at least for transitions between not too far away levels). We neglect virtual 
photons, whose contributions are suppressed by a. In QCD, it is also crucial 
to establish the size of these scales with respect to the typical hadronic scale 
Aqcd- By definition of heavy quark: m » Aqcd; however, mv » Aqcd 
holds possibly only for the lowest quarkonium resonances, while mv <~ 
Aqcd holds for the others. In the first case, called weakly coupled, the 
scale mv may be treated in perturbation theory (in this case v ~ a s ), in 
the last one, called strongly coupled, at the scale mv one cannot rely on an 
expansion in a s and non-perturbative techniques have to be used. 

In order to construct an EFT suitable to describe heavy quarkonium 
transitions, we do not need to resolve scales larger than or of the same or- 
der as the soft scale mv. These scales may be integrated out from the EFT. 
The relevant degrees of freedom differ in the weakly and in the strongly 
coupled case. If mv 3> Aqcd these are singlet, S, and octet, O, quarko- 
nium fields, gluons of energy and momentum mv 2 or Aqcd (sometimes 
also called ultrasoft) and photons. If mv ~ Aqcd, at scales lower than mv 
colour confinement sets in and the relevant degrees of freedom are singlet 
quarkonium fields (all gluonic excitations between heavy quarks develop a 
mass gap of order Aqcd with respect to the lowest state and are integrated 
out) and photons. Fields scale in accordance to the lowest, still dynamical, 
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energy scales, denning in this way the power counting of the EFT. 
The EFT Lagrangian has the following form: 
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where £ 7 is the part of the Lagrangian that describes the interaction with 
the electromagnetic field. In the case of magnetic-dipole transitions, the 
relevant terms in L 1 up to relative order v 2 are: 
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where e is the electric charge of the heavy quark. The Wilson coefficients V s , 
V , Va, Vb, V4 111 , Vi are in general functions of r and have to be determined 
by matching amplitudes calculated in the EFT with amplitudes calculated 
in QCD. The matching may be performed order by order in a s in the weak- 
coupling case, but has to be performed non-perturbatively in the strong- 
coupling one. In both situations, the matching can be done in two steps: the 
first one consists of integrating out hard modes of energy of the order of the 
heavy quark mass, the second one of integrating out soft modes of energy 
of the order of mv. The Wilson coefficients of the EFT will then have the 
factorized form: (hard) x (soft). For determinations of V s , V Q , which are the 
singlet and octet potentials, and Va, Vb we refer to 20 and references therein. 
The coefficients V| m and Vi have been studied in 21 . In the following, we 
summarize the main conclusions. 
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The coefficient V\ may be interpreted as the heavy quarkonium mag- 
netic moment. Its hard part is simply the sum of the heavy quark mag- 
netic moments: 1 + 2a s /(37r) + . . . . Moreover, it can be shown that, in the 
SU(3)ft a vour limit, to all orders in the strong-coupling constant, V\ does not 
get soft contributions. The argument is somewhat similar to the factoriza- 
tion of the QCD corrections in b — > ue~D e , which fixes the coefficient in 
the effective Lagrangian C e s = — 4Gf/\/2 V u b &Llti^L Ul^^l to be 1 to 
all orders in a s (see, for example, 22 ). This leads to the conclusion that the 
heavy quarkonium magnetic moment, if defined as the Wilson coefficient 
Vi, is just the sum of the heavy quark magnetic moments. Therefore, its 
anomalous part is small and positive and does not get any large low-energy 
and, in particular, non-perturbative contribution. This is consistent with 
a recent lattice determination 23 . Note that a large negative heavy quarko- 
nium anomalous magnetic moment has been often advocated in potential 
models to accommodate the results with the data. 

Reparametrization and Poincare invariance 24 ' 25 protect the coefficients 
V2 and V3. To all orders in perturbation theory and non-perturbatively it 
holds that V 2 = r 2 vj°^' /2 (vj°^ is the singlet static potential) and V 3 = 0. 
Since an effective scalar interaction would contribute to V3, we conclude 
that such an interaction is not dynamically generated in QCD for magnetic- 
dipole transitions. Again, a scalar interaction is often employed in potential 
models. 

The coefficient V% m is 1 up to possible corrections of order a 2 . In gen- 
eral, the coefficients of the operators of order 1/m 3 are not protected by any 
symmetry. V4 is 1 plus 0(a s ) corrections. In the weak-coupling case, this is 
the only 1/m 3 operator needed to describe magnetic-dipolc transitions at 
relative order v 2 . In the strong coupling regime, more terms are, in principle, 
necessary. Since a s (l/r) <~ 1 is no longer a suppression factor, more ampli- 
tudes shall contribute to the matching. These amplitudes will be encoded in 
the matching coefficients of the EFT in the form of static Wilson loop am- 
plitudes with field strength insertions of the same kind as those that appear 
in the QCD potential at order 1/m 2 26 . Also, they may induce new oper- 
ators in the EFT. A non-perturbative derivation of the EFT Lagrangian 
coupled to the electromagnetic field at order 1 /m 3 has not been worked out 
yet; note that to be of phenomenological impact such a calculation needs to 
be supplemented by lattice calculations of the relevant Wilson loop ampli- 
tudes. Since most of the potential model calculations rely on Lagrangians 
of the type of Eq. (2), we emphasize that, once cleaned up of the scalar 
interaction, they are suitable to describe the relativistic corrections for ra- 
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diative transitions between weakly-coupled quarkonia (mv 2 > Aqcd) only. 
This is a rather severe constraint since most of the radiative transitions 
(e.g. magnetic dipole hindered, electric dipole, ...) involve higher excited 
quarkonium states, which may be difficult to accommodate as Coulombic 
bound states. 

At present, allowed magnetic-dipole transitions between the quarkonium 
lowest states are those most reliably described in QCD. Moreover, it has 
been shown in 21 that for these processes, non-perturbativc contributions 
mediated by quarkonium octet states and ultrasoft gluons cancel at relative 
order v 2 . Therefore, at relative order v 2 , transitions like J/tjj —* r\d and 
T (15) — > r\b 7 are completely accessible in perturbation theory. An explicit 
calculation gives: 

r(J/V>^7 ?c7 ) = (1.5±1.0)keV, (3) 
T(T(15) Vb 7 ) = (gg^-y) (2-50 ± 0.25) eV. (4) 

The error in Eq. (3) accounts for the large uncertainties coming from 
higher-order corrections; we recall that corrections of order k^v 2 /m 2 af- 
fect the leading-order result by about 50%. Uncertainties may be reduced 
by higher-order calculations. The value given in Eq. (3) is consistent with 
the present experimental one . This means that assuming the ground-state 
charmonium to be a weakly-coupled bound state leads to relativistic cor- 
rections to the transition width of the right sign and size. The decay width 
r(T(15) — ► rjbj) depends on the r/b mass, which is unknown. The quoted 
value corresponds to a rjb mass of about 9421 MeV. 

We conclude with few possible developments of the work discussed here. 
For magnetic-dipole transitions between higher resonances, the completion 
of the non-perturbative matching of the relevant operators in the EFT at 
order 1/m 3 will be needed, possibly integrated by lattice calculations of the 
relevant Wilson-loop amplitudes. For magnetic-dipole hindered transitions 
of the type T(35) — > f]b7, since the momentum of the emitted photon 
is comparable with the typical momentum transfer in the bound state, 
one cannot rely on the multipole expansion of the external electromagnetic 
field. In this case, one may, in principle, exploit the fact that the typical 
momentum transfer inside the r\b is much larger than that one inside the 
T(35). A suitable treatment has not been developed yet. For electric-dipole 
transitions, much of the study still remains to be done. In the weak-coupling 
regime, octet contributions may not vanish; they can be worked out as in 
the case of the magnetic-dipole transitions. However, most of the electric- 
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dipole transitions may need to be treated in a strong-coupling framework. 
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